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The task of the research was to build and analyze a model of the dynamics of 
pollution of a flowing reservoir and systems of communicating reservoirs 
with and without taking into account water filtration in the soil as a result of 
external sources (effluents from industrial enterprises). This work was aimed 
at studying the change in temporal dynamics, taking into account the 
concentration of impurities in the volumes of three reservoirs during the 
periods of discharge and completion, lasting 30 days. Numerical 
experiments were carried out for various flow rates and compositions of 
filtration coefficients to study the relaxation times of pollution in the system 
of reservoirs. Also, software was developed that analyzes the change in the 
concentration of impurities in the system of reservoirs. As a result, it was 
found that the pollution pattern is more dependent on the topology of the 
watercourse network. 
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1. INTRODUCTION 


Many scientific research papers consider mathematical modeling of water systems [1]-[5]. The 
analysis of the considered literary sources showed the possibility of conditionally dividing the main 
mathematical models of water systems into 2 classes: those used in forecasting and those used for 
monitoring. Mathematical models of water systems have a number of features. First of all, this is due to the 
very object being modeled-water. Models describing the water environment are multi-criteria since they 
provide both economic and environmental needs of society. When setting such tasks, there is a need for an 
integrated and versatile approach to the study of the problem [6], [7]. Another feature of the water system is 
its dynamism, i.e., any process in the environment should be considered in time. Research papers [8]-[10] 
show that mathematical models used to study the water environment are complex and diverse in their 
structure. To build models of this type, it is very important to correctly formulate the problem statement and 
generalize the available initial data, then choose the optimal solution method. 
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Analysis of literary sources shows that, despite a large number of studies, the known models 
describing the dynamic distribution of pollutants in the water system are complex and limited. Theoretical 
concepts of impurity spreading ways in the water environment and distribution in the system of water 
reservoirs require further development and improvement [11]-[13]. Today, the study of the degree of 
pollution in reservoir systems has received great practical interest. The results of such studies are widely used 
to solve various problems in ecology [14]-[18]. This approach involves the use of simplified samples for the 
analysis of dynamic pollution of running water reservoirs in neighboring systems of reservoirs, considering 
options with and without taking into account water filtration in the soil under the influence of external sources 
(discharges from industrial enterprises). The next step for the development of the model is the application of an 
analysis of the change in the concentration of impurities in the system of co-correlating reservoirs under 
consideration, taking into account the diffusion process of impurities in water bodies and soil filtration, i.e., in 
the dynamic model, typical periods of impurity propagation within each reservoir should be provided. 


2. RESEARCH METHOD 

Industrial emissions can be distributed in a system of interconnected storage facilities, the result of 
which for the most part is a rather complex pollution situation [19]-[21]. This article examines a system 
consisting of three reservoirs. In the coastal zone of one of the reservoirs, there izs a source of pollution with 
a given intensity Figure 1. 

At the beginning, the condition of total expense balance by flows in the system was offered. This 
condition means, in fact, neglecting the water filtration system in the soil. Further, in order for the 
development of the model to be directed to the dynamic analysis of the concentration of additives in the 
system of communicating reservoirs, taking into account diffuse impurities in reservoirs and filtration in the 
soil, it was necessary to introduce the characteristic times of impurity propagation in the volumes of each 
reservoir in the dynamic model [21]. Figure 2 shows the flow diagram in the system of three communicating 
reservoirs, taking into account filtration. 
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Figure 1. The flows scheme in the system of three Figure 2. Flows scheme of the three connected water 
co-correlating reservoirs without taking into account reservoirs system with the filtration consideration 


water filtration in the soil 


The article highlights the main points that are associated with the dependences of the following 
characteristics: on the relaxation time in each of the reservoirs, depending on the duration of the discharge, as 
well as the coefficients of water filtration in the soil [21]. Therefore, the main goal of our research was to 
determine the dependence between the concentrations of impurities in each of the reservoirs on the intensity 
of filtration in the channels connecting these reservoirs in a given period of time [22]. 


2.1. A case of no filtration in the soil 
Let’s to write down the system of differential equations for calculating the balance of pollution of 
the reservoir with stable impurities, provided that the condition of filtration in the soil can be neglected, 


V, = 1(t)-s*q 


dcz 


V2- = Ss) * c — (s = s1) * c2 (1) 
d 
T= s * G + (S— 51) * C2 — S * C5 
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where c4, C2, C3 are mean impurity concentrations relatively in the 1*, 24 and 3" reservoirs, defined in kg/m’; 
S, S4, S2, S3 are water flows, defined in m*/sec; Vi, V2, V3 are water volumes, defined in m; I is discharges 
intensity, defined in kg/sec; f is the time, defined in sec. The following conditions conforming to a total flow 
balance are accepted for the system. 


s1 + s2 = s; s2 = s3; s1 + s3 =s, (2) 
Let us consider the impurity concentration dynamics in each of the water reservoirs in the case of 


discharge impact with constant intensity I(t) = I = constduring some period T. In this case, it is possible to 
present an approximation function of the impact discharge in the form of a rectangle as shown in Figure 3. 


` 


T 


Figure 3. Scheme of the peak discharge 


The impurities concentrations in the water reservoirs are c,(0) = C1(9), C2(0) = C2(9), C3 (0) = C3(0) - in the 
case of the initial discharge moment. 


2.2. A case of the filtration availability through the soil 
It is possible to transform a dynamic model (1) of pollution of reservoirs with impurities that have a 
stable composition and take into account the filtration process: 


V, = = I(t) — So #0, = I(t) — (S1 + s2) * c (I) 
d 
Vp 2 = Q, * S1 * C1 — Q; * S1 * C2 UD (3) 
d 
Va = Py * S2 * Cy + P3 * P, * S1 * C2 — (P, * Pz * S1 + S1S1C1 + Y,S2) * C3 (I) 


where C1, C2, C3 can be considered as the impurity concentrations relatively in the 1%, 2" and 3" reservoirs, 
kg/m*; So, S1, S2, Sg can be considered as the water flows, defined in m*/sec; Vi, V2, V3 are water volumes, 
defined in m*; I is discharges intensity, defined in kg/sec; tis the time, defined in sec; @,,@2,3 are the 
coefficients of filtration. Conditions (4) were accepted for the system (3). 


S1 + S2 = So; Sk = $1 * hz * S1 + hz * 82 (4) 


In the initial discharge moment, the impurities concentrations in the water reservoirs were, 
respectively: 


c& (0) = C1(0)» C2 (0) = C2(0)) C3 (0) = C3(0) (5) 


in the article [23], it was was considered a solution of the system (3). The next section gives the results of 
numerical experiments for both cases. 


3. NUMERICAL EXPERIMENT RESULTS 
Figures 4 and 5 show graphs that fix the temporal dynamics of the concentration of impurities 
during and after the discharge, neglecting water filtration in the soil. It can be seen here that the average 
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values of the concentration of impurities in individual tanks during the discharge process begin to increase, 
and then, having reached certain values, they reach a plateau. Consequently, a constant total concentration 
reaches equilibrium in each of the reservoirs, provided that the discharge has a constant intensity for a long 
time. At the same time, the dynamics of concentrations in individual water bodies can have significant 
differences during a short discharge duration. Further, we can notice that after the start of the discharge, the 
impurity concentration in the first reservoir increases much faster compared to the same concentrations 
c,(0) = 0, c2(0) = 0, c3(0) = 0 in the remaining reservoirs. For example, the following time intervals 
were considered: 4,65-10° sec, 4:10° sec and 1-105 sec. Consequently, the experimental data at T=4,65-10° 


impurity concentrations in water bodies when discharged into all three reservoirs were 99,991 kg/m?, 99,653 
kg/m?, 99,413 kg/m’, respectively. 
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Figure 4. Graphs of changes in time of impurity concentration in 3 reservoirs, occurring during and after 
discharge at T=4,65-10 5 sec 
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Figure 5. Graphs of changes in time of impurity concentration in 3 reservoirs, occurring during and after 
discharge at T=4-10 5 sec 


At the same time, impurity concentration curves (C,,C,,C3) were calculated for a certain set of 
filtration coefficients in all reservoirs. The following are examples of filtration coefficients: 
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1. Ø, =0,2, p, =0,6, Q; =08; 2. P, =1, P, =1, p, =1; 


3. Ø, = 0,2, 
P, =0,2, p, =0,2; 4. p, =0,8, 


P, =0.8, p, =0,8; 5. Ø, =0,2, p, =0,8, P, =0,6; 6. p, =0,8, p, =0,6, p, =0,2. 


Discharge duration options from 200000 to 800000 sec were considered. Below you can see some 
results of the numerical experiment. For example, Figures 6 and 7 show graphs of the time dependence of the 
concentration in the |-st reservoir at @; =0,2, d2 =0,8, @, =0,6 for various T1, T2, T3, T4. 
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Figure 6. “C — t” curves of the impurity dynamics in the first reservoir in various discharge durations T1, T 2, 
T3, T4 (hı =0,2, 2 =0,8, p3 = P3 = 0,6) 


Tl T2 T3 T4 T 
Figure 7. Graph of the time dependence of relaxation on the duration of the discharge in the first reservoir at 


$ı =0,2, $2 = b2 =0,8, 3 = p, =0,6 


Next, relaxation periods T were determined for the duration of each discharge, characterized by 
the fact that the concentration during this period decreases sharply after the completion of the discharge from 
Cmaxto 0.05Cmax. The experimental data obtained in the course of numerical experiments showed that the 
relaxation periods in the 1l-st reservoir are similar and have a value of about 150000 sec, taking into account a 
certain discharge amplitude. With an increase in the duration of the discharge, the relaxation period increased 
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in the second reservoir, which was considered with similar filtration coefficients and discharge duration 
values calculated in the first reservoir. Summarizing the data obtained from the graphs, we can say that the 
dynamics of the flow duration in the first reservoir does not depend on filtration, just as the relaxation phases 
in any filtration installations do not have any dependence on the flow duration in the first two reservoirs. But 
at the same time, the third reservoir has a relaxation phase, which directly depends on the given filtration 
coefficients, as well as on the duration of the discharge. 


4. APPROVAL OF SOFTWARE SYSTEM FOR ANALYSIS OF TOXIC IMPURITIES 
CONCENTRATIONS IN THE SYSTEM OF COMMUNICATING WATER RESERVOIRS 

Modern research-practical environmental problems put qualitatively new requirements on the 
computer processing of information. Similar problems are characterized by labor-intensive processing of 
input and output information, in this connection, the role of computational simulation and experiment 
increases [24], [25]. This section considers questions of creation and approval of software systems for 
analysis of changes in the impurity concentration in the systems of communicating water reservoirs. Earlier 
we considered systems consisting of three running reservoirs with the pollution source of specified intensity 
on the bank of one of the reservoirs [25]. Mathematical models that describe the dynamic pollution of 
reservoirs with various stable impurities, both with and without filtration in the soil, are described in these 
works. Decisions, worked out in the course of the research objective, were directly used during development 
of the software system in DELPHI environment. The software system is assigned for automation of 
processing and computation using specified initial data. The program has means of dialogue interaction with 
a user at fulfilment of special actions, and convenient interface. A project, implemented in DELPHI 
environment, consists of eight interrelated forms. After the program start, a window with information about 
software system appears as shown in Figure 8. This window is the main one, and it is used for transit to other 
windows of our application. After choice of one of possible options in the list, the required window of the 
program appears. Figures 9 and 10 show windows of programs for two studied cases. 
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Figure 8. Basic form 
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Figure 9. Window for analysis of the impurity concentration change without the filtration 
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Figure 10. Window for analysis of the impurity concentration change with filtration 


Each of the above windows contains 4 buttons: “compute”, “mathematical model”, “diagram”, and 
“close”. For operation of the program, in the left field of the windows, it is necessary to set data, defining the 
duration of the discharge, concentration of pollutants in all water reservoirs in the initial moment of time, 
capacity of the reservoirs, rate of water by flows, and intensity of the discharge, moreover, values of the 
filtration coefficients are additionally put in for the window of Figure 8. The output data of the program are 
mean concentrations of the impurity in each of the reservoirs in the specified time, presented in the right field 
of Figures 7 and 8, and also diagrams of temporal dynamic concentration of the impurity in the water 
reservoirs both during and after the discharge. Figures 11 and 12 show diagrams of temporal dynamic 


concentration of the impurity in three water reservoirs during and after the discharge at T=4,65-10 5 sec with 
consideration and without consideration of the filtration. 
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Figure 11. Diagram of temporal dynamic concentration of the impurity in three water reservoirs during and 


after the discharge at T=4,65:10 5 sec without consideration of the filtration 
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Figure 12. Diagram of temporal dynamic concentration of the impurity in three water reservoirs during and 


after the discharge at T=4,65-10 5 sec with consideration of the filtration 


Obtained results allow easily to determine the values of the impurity concentrations in each of the 
reservoirs at any moment during and after the discharge. “Mathematical model” button in the windows 
(Figures 9 and 10) gives a graphic theoretical presentation of the studied research objects. For convenient 
review of the obtained results, there are auxiliary tools in the program. Computing results are automatically 
stored in files of “Notes” in the standard program. The graphic part of the program is implemented by 
Diagram component (TChart). This strong and high-powered component was developed by David Berneda, it 
allows construct complex diagrams on the basis of different data. The operability of the program is confirmed 
by the fact that the experiments were carried out with given initial conditions obtained from analytical 
solutions, which, in turn, were correlated with the numerical data of the calculated results using application 
programs such as Mathcad Prime 2.0, and MATLAB. Analysis of the software system allows to make the 
following conclusions: 

— The program allows operating by large bodies of input data; 

— Results, obtained during operation with the program, differ by high accuracy; 

— The software system has sufficient clear visualization, as it was performed in the environment of visual 
programming DELPHI; 

— The program has high performance and operability; 

— lso, the program can be used in the process of modeling tasks that solve the problems of the distribution 
of industrial discharges in the system of communicating reservoirs. 


5. CONCLUSION 

The research analysis allows to conclude that the picture of the propagation of pollutants contained 
in waste waters entering water systems from production depends on the type of industrial enterprise. 
Accordingly, when studying the interaction of impurities in the water environment, it is necessary to take into 
account their qualitative composition. The simplified model for analyzing the relaxation propagation of the 
impurity in the system of communicating water reservoirs has been submitted. Approval of the model was 
carried out for the systems, which consist of three running reservoirs both with consideration and without 
consideration of the water filtration in the soil. The system of equations for balance of the water reservoirs 
pollution by the stable impurity was derived and solved. The temporal dynamic concentration of the impurity 
in three water reservoirs during and after the discharge, during certain period, in the systems without 
filtration through the water reservoir bottom was studied. The similar model for dynamic pollution of the 
systems of communicating reservoirs with consideration of the water filtration in the soil was constructed. 
The analysis of dynamic concentration of the impurity in the system of communicating reservoirs with 
consideration of the impurity diffusion in the reservoirs and filtration in the soil was also carried out. 
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Water systems function and develop in time and space under the influence of internal laws. The 
existing change tendencies in the system parameters remain in the form of a trace over a certain period. With 
this approach, it becomes possible to identify not only the qualitatively correct tendency of the system, but 
also to obtain calculated dependencies for a reliable forecast. Relaxation models are highly flexible and 
adaptable. The numerical experiments at different discharge durations and compositions of the filtration 
coefficients were carried out to study the pollution relaxation times in the system of water reservoirs. It was 
revealed that the relaxation time in the first and second reservoirs does not depend on the discharge duration 
at any filtration compositions, however, this dependence is observed for the third reservoir. The pollution 
relaxation time dependence in the third reservoir is explained by location of this object in the system of 
communicating water reservoirs. It was established that the pollution picture to a large extent depends on the 
scheme of the water flows net. The pollutant particles’ propagation in the water environment depends on 
their size and chemical nature; therefore, relaxation occurs at different rates, being a multi-step process. 
Various systems are characterized by their own relaxation features due to the nature of the interaction 
between the particles of the system. 

During modeling of the industrial discharge propagation in the system of communicating water 
reservoirs, the questions of creation and approval of the software system to analyze the impurity 
concentration change in the system of these water reservoirs were considered. Results obtained during 
operation with the program can be applied in researches of the industrial discharges propagation modeling in 
the system of communicating water reservoirs. It is possible to use our models to calculate these effects, 
while adequately predicting the change in the concentration of an admixture in various systems of connected 
reservoirs. Also, based on this approach, it is possible to calculate the maximum values of concentrations for 
forecasting the time, the duration of which is necessary when determining the excess of the maximum 
permissible concentration after the discharge of impurities into reservoirs. 

The results of the calculations led to the fact that there is some time delay in the dependence of 
concentration on time for the impurity concentration in each of the reservoirs. However, after the discharge is 
completed, there may be a temporary increase in the concentration of impurities in the event of a delay in a 
water body remote from the point of direct discharge, even in the case of a small amount of discharge. 
Observations of temporal changes in pollution in the case of discharge and relaxation confirmed numerical 
experiments, which showed the adaptability of the developed model. This suggests that this method can be 
recommended for research, analysis, and forecasting of dynamic environmental pollution in industrial areas. 
Thus, the use of mathematical methods, models and special software products allows to study the features of a 
water body and its sources of pollution, carry out extensive monitoring and make a qualitative comparison with 
indicators of previous years, and assess the functioning and development of the water system in the near future. 
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